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ABSTRACT 

Using 2D and 3D simulation, we study the "robustness" of the double detonation scenario for Type 
la supernovae, in which a detonation in the helium shell of a carbon-oxygen white dwarf induces a 
secondary detonation in the underlying core. We find that a helium detonation cannot easily descend 
into the core unless it commences (artificially) well above the hottest layer calculated for the helium 
shell in current presupernova models. Compressional waves induced by the sliding helium detonation, 
however, robustly generate hot spots which trigger a detonation in the core. Our simulations show 
that this is true even for non-axisymmetric initial conditions. If the helium is ignited at multiple 
points, the internal waves can pass through one another or be reflected, but this added complexity 
does not defeat the generation of the hot spot. The ignition of very low-mass helium shells depends 
on whether a thermonuclear runaway can simultaneously commence in a sufficiently large region. 
Subject headings: hydrodynamics - nuclear reactions, nucleosynthesis, abundances - shock waves - 
supernovae: general - white dwarfs 



1. INTRODUCTION 

In the double detonation scenario, a carbon-oxygen 
(CO) white dwarf accumulates a layer of helium through 
accretion from a helium main sequence star. The accre- 
tion rate, typically on the order of 1O~®-1O~^M0 yr~^ 
(e.g., Taam 1980; Woosley & Weaver 1994), and the mass 
of the white dwarf determine the critical thickness of the 
helium layer. Generally, a low accretion rate allows more 
helium to accumulate before a thermonuclear runaway 
commences near the base of the helium shell, the tem- 
perature of which is determined by gravitational com- 
pression, convection and eventually nuclear burning. It 
is often assumed that a detonation in the helium induces 
a secondary detonation in the CO core, thus producing 
a Type la supernova. 

If and how a secondary core detonation develops is 
ambiguous, however. One possibility is that the core is 
ignited directly at the core-shell interface when hit by 
the helium detonation front. This is sometimes called 
"direct drive" , "edge-lit scenario" or "prompt detona- 
tion" . It has long been realized that this does not always 
work and that its realization depends on the altitude at 
which the helium detonation starts (e.g., Livne & Glas- 
ncr 1990; Bcnz 1997; Garcia-Senz et al. 1999). A helium 
detonation directly at the interface almost certainly fails 
to start a detonation in the core, but evolution models 
of accreting white dwarfs suggest that the most proba- 
ble place for a detonation to start could be higher up 
(Woosley & Kascn 2011). One of the aims of the present 
work is to determine whether a helium detonation ignited 
at a point in the hottest helium layer would lead to core 
detonation by direct drive. 

If direct drive fails, the core may still ignite, with con- 
siderable delay, after the detonation wave has consumed 
the entire helium layer, through compressional waves 
which focus inside the core. In the simplest (and most 
unrealistic) scenario, the star is spherically symmetric 
and the helium is ignited instantaneously in a spherical 
shell. This induces a radially inward propagating com- 
pressional wave that converges at the exact center of the 



star, where, depending upon the resolution of the study, 
it always ignites the core (e.g., Livne 1990). 

Since convective mixing is not able to keep the temper- 
ature at a given radius constant before the runaway, it 
is more likely that helium ignition starts in one or more 
isolated regions with sizes of perhaps a pressure scale 
height (Woosley & Kasen 2011). If the detonation starts 
at a single point, the sliding helium detonation creates 
an oblique (with respect to the core-shell interface) shock 
wave that converges inside the core at the axis of sym- 
metry and potentially triggers an off-center core detona- 
tion (Dgani & Livio 1990; Livne & Glasncr 1990, 1991; 
Livne & Arnett 1995). In their simulations of double- 
detonation supernovae. Fink et al. (2007) considered a 
variety of ID and 2D helium ignition scenarios, confirm- 
ing the high ignition potential of the inward propagating 
shock wave. They showed that core detonation result- 
ing from the detonation of the helium in a single spot is 
very likely even for low mass helium shells, for instance 
O.OO35M0 in combination with a 1.385OM0 core (Fink 
et al. 2010), and for core masses as low as O.45M0 (Sim 
et al. 2012). 

In a scenario with less or no symmetry, if helium det- 
onations commence at multiple points, for example, the 
secondary detonation of the core is less certain. In their 
SPH simulations, Garci'a-Senz et al. (1999) first explored, 
using coarse resolution, the possibility of starting deto- 
nations at multiple points in the helium layer. In partic- 
ular, they studied a case with asynchronous ignition at 
five points, finding that core detonation could still result 
from the collision of the various helium detonations. 

If the core does not ignite, the helium explosion alone 
could produce a sub-luminous, "point la" , supernova 
(Bildsten et al. 2007). This raises the question as to 
how robust the triggering of the secondary core detona- 
tion is. If it is robust, a detonation of the helium shell 
alone might be very rare or even impossible. 

We here revisit the double detonation scenario using 
some of the white dwarf models generated with KEPLER 
and presented in Woosley & Kascn (2011). Following a 
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description of the numerical methods in §2, we first dis- 
cuss the results of simulations that test the viability of 
direct drive in §3.1. We then (§3.2) study sliding helium 
detonations that lead to a hot spot inside the core, with 
the helium detonation started at one, two or three points. 
The initiation of detonations in lightweight helium shells 
is investigated in §3.2.6. In §3.3, we present nucleosyn- 
thesis yields from runs following the complete detonation 
of white dwarfs. We end with a discussion and summary 
in §4. 

2. NUMERICAL METHODS 
2.1. Hydrodynamics 

For our multi-dimensional studies we use the Eule- 
rian hydrodynamics code CASTRO (Almgren et al. 2010; 
Zhang et al. 2011) to solve the equations for compressible 
fluid dynamics in combination with a nuclear reaction 
network. We employ an equation of state based on the 
Helmholtz free energy (Timmes & Arnett 1999; Timmes 
& Swesty 2000) and a 19 isotope network to compute 
energy generation and nucleosynthesis (the same as in 
KEPLER, see §2.2 for details). In many simulations, we 
make use of CASTRO 's capabilities for adaptive mesh 
refinement (AMR). The initial models for the multi-D 
CASTRO simulations presented here were constructed 
with the implicit Lagrangian hydrodynamics code KE- 
PLER in ID (Weaver & Woosley 1978; Woosley et al. 
2002). The evolution of several helium-accreting white 
dwarfs was followed with a model for convection based on 
mixing-length theory, see Woosley & Kascn (2011). The 
CASTRO simulations start from KEPLER snapshots 
taken before a thermonuclear runaway sets in. Their key 
properties are listed in Table 1. 

In Eulerian simulations, adding a relatively small 
amount of hot gas to a cold zone can raise the average 
temperature above the threshold for ignition, thus ignit- 
ing the whole zone. Without any correction, the speed 
with which a "detonation front" propagates would thus 
be dependent on the size of the zones and driven entirely 
by artificial advection. To simulate detonation more re- 
alistically, we introduce a delay timer that suppresses 
the burning in a particular zone for half a sound crossing 
time. If, at any time, the necessary conditions for nuclear 
burning are not satisfied, the timer is set to zero. If they 
are satisfied, the timer starts counting and when half a 
sound crossing time (Courant time) has passed, nuclear 
burning processes are switched on in this particular zone. 

In simulations of the whole star, the density of the ar- 
tificial external medium is set to 1 g cm~'^ and given the 
same temperature as the outer edge of the helium shell. 
This means that as the shock of the helium detonation 
breaks out of the star, it is accelerated by a very steep 
pressure gradient. As the time step of our simulations is 
limited by the largest fluid velocities or sound speeds in 
the entire domain, the calculations become very expen- 
sive. To contain the detrimental effect on the time step, 
we apply a velocity cap of 1.5 x 10^ cm s~^ in the ambient 
medium. To distinguish the ambient from stellar mate- 
rial at all times of the simulation, we give it a unique 
composition of 100% ^^O. We used no velocity cap and 
lower ambient densities (starting with 10^^ g cm~^ and 
ending with 10~^ gcm~^ in the final computational box) 
in runs where the expansion of the ejecta was followed 
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Fig. 1. — Variation of the final composition for a peak temper- 
ature of 3 X 10^ K when the density in the simulated explosion is 
varied from 10^ to lO^gcm""^. The solid and dashed lines rep- 
resent the results obtained with the full network and the table, 
respectively. 
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Fig. 2. — Variation of the final composition for a constant den- 
sity of lO^gcm"^ when the peak temperature in the simulated 
explosion is varied from 1 to 4 X 10^ K. The solid and dashed lines 
represent the results obtained with the full network and the table, 
respectively. 

(§3.3). 

We generally make use of symmetries to reduce the 
computational load. With only one detonator, the prob- 
lem is axisymmetric and its solution can be computed 
in 2D. For a simulation with two synchronous spherical 
detonators it is sufficient to do a 3D calculation of only a 
quarter of the star, with mirror symmetry at two bound- 
aries (the a; = and y = boundaries in our simula- 
tions). Two asynchronous detonators require a half-star 
simulation with mirror symmetry at one boundary (the 
X ~ boundary in our simulations). 

2.2. Nuclear Energy Generation and Nucleosynthesis 

Energy generation and approximate nucleosynthesis 
were calculated in all cases using a 19 isotope network 
with constituent species as defined by Weaver & Woosley 
(1978) and nuclear reaction rates updated to current 
values (Woosley & Heger 2007). Screening corrections 
were implemented for all rates. This network includes an 
"alpha-chain" {Z = N = 2n, n = an integer) of nuclei 
from ^^C through ^^Ni, plus the species n, p, '^He, ^^N, 
and ^"^Fe. Simulated steady state links also approximate 
the presence of the odd Z isotopes from ^'^Na through 
^^Co (i.e., isotopes with Z = 2n+l, N = 2n+2). For this 
study, the photodisintegration of ^^Ni to ^^Fe -I- 2p was 
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TABLE 1 
White dwarf models 



Model'^ 


-^core 




''^corc 




'"('j^max ) 


^max 


p(r(Tmax)) 


P(rcorc) 


p{r = 0) 




[Mq] 


[Mq] 


[km] 


[km] 


[km] 


[10« K] 


[10*5 gcm-^] 


[106gcm-3] 


[lO^gcm-3] 


A (SB) 


0.801 


0.143 


4081 


2080 


4267 


2.06 


1.19 


1.66 


2.21 


B (lOB) 


1.000 


0.082 


3616 


1370 


3765 


2.42 


1.30 


1.93 


5.67 


C (IOC) 


1.000 


0.091 


3518 


1380 


3624 


2.31 


1.63 


2.18 


6.02 


D (lOHC) 


1.002 


0.045 


4161 


1260 


4187 


2.80 


0.717 


0.932 


4.24 


E (lOHD) 


1.001 


0.078 


3696 


1350 


3724 


2.36 


1.50 


1.77 


5.33 



^ The model identifiers in Wooslcy & Kasen (2011) are written in parentheses 

^ Approximate helium shell thickness, with the outer boundary defined where the density is 1/100 of the value 
at the core-shell interface 



not allowed which meant that electron capture and ^''Fe 
were essentially omitted from the network. Furthermore 
^H, ■^He, ^*N, n, and p never achieved any significant 
abundances, so the effective network really just contained 
the 13 alpha-chain isotopes. Past studies (Timmes et al. 
2000) have shown that such a limited network gives good 
agreement with the results of much larger networks pro- 
vided one is interested only in bulk nucleosynthesis and 
energy generation. 

The network itself was employed in all carbon-rich 
zones always and, in most cases, for the helium shell. 
In some cases, however, it was computationally expedi- 
ent to use a table that was prepared off-line using the 
network. Those cases included studies of helium shell 
detonation where the whole star was carried (§3.2), up 
until the generation of a hot spot in the core. The sim- 
ulations for which detailed yields are given (§3.3) make 
use of the full network from beginning to end, as do the 
local simulations of direct drive (§3.1). 

The use of a table here was novel and warrants some 
discussion. Unlike isobaric burning in e.g., carbon defla- 
gration supernovae, the table here has to work reliably 
in dynamic situations where both the temperature and 
density are rapidly varying and where the outcome is 
sensitive to the time history of both. Using the network, 
a large table was generated using the results of many off- 
line studies of helium burning at constant temperature 
and density. Values of temperature in the range lO*'^ - 
109-65 (AlogT = 0.05) and density in the range lO^^ - 
10''-^gcm~^ (Alogp = 0.10) were included. Results of 
the burning - the composition, average binding energy 
per nucleon, and time derivative of the helium mass frac- 
tion, were sampled for 50 values of helium mass fraction 
ranging from 0.02 to 1.0 in steps of 0.02. To preserve 
accuracy, the rate equation for the helium mass fraction 
was divided by the leading order term in the 3a reaction 
rate, i.e., 



dt 



tabic 



\ dt 



rp3 

-^9 



calc 



p^Xl exp 



-4.403 



(1) 

This resulted in a slowly varying function of Tg, p, and 
Xa in the table. 

When using the table in a multi-dimensional simula- 
tion, triple variable interpolation in logp, logTg, and Xnc 
gave the current composition and dXijo/dt. The current 
time step and dXn^/dt gave a new helium abundance at 
the end of the step. Another call to the table using the 
new helium abundance gave a revised composition for 
the other species, and the change in composition gives a 



change in nuclear binding energy, hence an energy gen- 
eration, and the process continued. 

The shortcoming of this approach is that the entries in 
the table were generated at constant temperature. This 
neglects additional variables because a given tempera- 
ture, density, and helium abundance can be reached by 
a variety of histories and the instantaneous abundances 
of the other nuclei are sensitive to which history was 
followed. Some additional rules are imposed to avoid 
unphysical and poorly behaved solutions. The helium 
abundance was only allowed to decrease, i.e., dXnc/dt 
had to be negative, or no composition change or energy 
generation is allowed. Further, no modification of the 
composition occurred unless both the abundance of ^^Ni 
and the nuclear binding energy increased during a time 
step. When these conditions were met, a change in com- 
position was allowed, but the new composition was not 
just taken from the table. Instead a fraction of the new 
composition at the given temperature and density was 
added to the old one with that fraction determined by 
the change of the helium abundance during the step. Es- 
sentially this procedure partitions the nuclear evolution 
for an arbitrarily varying temperature and density his- 
tory into a series of small steps taken at constant tem- 
perature and density. For the conditions of the problem, 
the requirement that the time step be a small fraction of 
the Courant time proved adequate to preserve accuracy. 

To demonstrate the validity of the procedure, a series 
of simulated explosive helium burning studies were car- 
ried out off line using both the network and the table. 
Starting with a composition of pure helium, matter was 
expanded from a peak temperature, Tgo at constant den- 
sity on a time scale t such that 



Tg(t) = Tgo exp(-t/r). 



(2) 



and the final composition examined when the tempera- 
ture had declined below 10* K. The results are shown 
in Figs. 1 and Fig. 2. Figure 2 shows the variation of 
the final composition for a constant density of 10^ gcm^'^ 
when the peak temperature in the simulated explosion is 
varied from 1 to 4 x 10^ K. In general the agreement be- 
tween the network results and those obtained using the 
table is quite good. There is a substantial deviation of 
up to 50% in the final helium abundance for the higher 
temperature explosions. Using the table gives a smaller 
final abundance for helium because burning at a constant 
low temperature gives a smaller mean nuclear mass and 
charge than explosive burning cooling from a higher tem- 
perature to the same helium mass fraction, density and 
temperature. Lighter nuclei have larger cross sections for 
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TABLE 2 

Direct drive: masses of the helium detonation 
spheres at core contact 



model 


^detonator 

[km] 


[km] 


-^^H c- s p licr c 

[lO^gcm^^] 


works /fails 


A 


4270 


189 


3.89 X 10^** 


fails 


A 


4470 


389 


2.80 X 10^^ 


fails 


A 


4670 


589 


8.28 X 10^9 


works 


B 


3765 


149 


2.17 X 10^** 


fails 


B 


3815 


199 


4.75 X 10^* 


fails 


B 


3865 


249 


8.71 X 10^® 


works 


C 


3624 


106 


9.40 X 10^^ 


fails 


C 


3636 


118 


1.26 X 10^8 


fails 


c 


3649 


131 


1.68 X 10^8 


works 


c 


3674 


156 


2.75 X 10^8 


works 


c 


3724 


206 


5.91 X 10^8 


works 



Distance from the core-shell interface to the detonator 

center 

^ Defined as twice the mass of the helium contained in 
the half-sphere of radius Ar below z = rdctonator 



a-capture. One could improve the fit by adding a fourth 
independent variable, the mean nuclear charge (exclud- 
ing helium itself), but the present approach was deemed 
adequate for this study. Some test runs at constant tem- 
perature and density using the table and the network 
gave near perfect agreement showing that the table had 
been correctly implemented. 

Figure 1 shows the results for a series of explosions with 
constant peak temperature (3 x 10^ K) calculated for a 
variety of densities. For ease of plotting and because it 
most directly affects the energy generation, the species 
Si, S, Ar, and Ca have been combined into "intermediate 
mass elements" (IME) and the elements Ti, Cr, Fe, and 
Ni have been combined into "iron group elements" (IGE) . 
^^Ni is also shown separately. The overall agreement 
between network and tables is again quite good. 

3. RESULTS 
3.1. Direct Drive 

Here we consider the possibility of a helium detona- 
tion wave that traverses the core-shell interface beneath 
its origin, directly igniting the core at the outer edge. 
Taking the (0.801 + O.143)M0 white dwarf Model A as 
initial condition, we did a series of 2D simulations that 
cover a local region including the core-shell interface. 
The domain range in most cases is < i?8 < 1 and 
3.3 < zg < 5.3 in cylindrical coordinates (i?, z), the core- 
shell interface being at the spherical radius rg = 4.08. 
The grid resolution in most cases is 0.977km (no AMR). 
The gravitational field is taken to be static and the 
boundary conditions are usually reflecting on all sides 
(as the detonation proceeds supersonically, the choice 
of boundary conditions is not critical). All simulations 
employed a full 19-isotope network which covers helium, 
carbon and oxygen burning. The detonation spots used 
to start the detonation had a radius of 50 km, a tem- 
perature of 2 X 10^ K at the center and 1.8 x 10^ K at 
the outer edge (linearly decreasing), a radial velocity of 
8000 km s""'^ at the outer edge (increasing linearly from 
zero) and twice the ambient density. 

Figure 3 depicts snapshots of several cases. In the stan- 
dard case (a), the helium detonator is centered at the ra- 
dius at which the temperature is maximal, 189 km above 



the core-shell interface. While the detonation creates an 
indentation in the interface, there is only little carbon 
burning and the detonation does not penetrate the core. 

With the helium detonation initiated 200 km higher, 
case (b), the detonation penetrates into the core. The 
front is unstable though, developing dents filled with 
compressed material that fails to ignite. While it is pos- 
sible that core ignition ensues from a detonation that 
fragments, we did not follow such pathological cases in 
our study and dismissed them as failures. In a similar 
simulation with half resolution (1.95km), case (c), the 
front is stable. The oxygen starts to burn 1130 km be- 
low the interface (at 72% of the core radius), where the 
density exceeds 5.3 x lO^gcm"'^. 

With the helium detonation initiated 400 km above the 
temperature maximum, case (d), the detonation pene- 
trates the core at our standard resolution, not showing 
any instabilities. We followed this detonation in a global 
simulation that encompasses half of the star, with the 
same resolution («1 km) as in the local simulation in a 
region that extends from the detonation spot to the cen- 
ter of the star {R% < 1 and zg < 5). This detonation 
stably proceeds all the way to the center. Oxygen burn- 
ing starts 830 km below the interface (at 80% of the core 
radius) . 

The fronts of detonations started at different altitude 
have different curvatures when they hit upon the core- 
shell interface. This raises the question as to whether 
the curvature of the detonation front plays a key role. A 
detonator in the form of an oblate spheroid fares better 
than a spherical one, compare panels (e) and (a). How- 
ever, the detonation front starts out considerably hotter 
than in the spherical case, which is possibly related to 
the details of how the detonation was started. 

We also considered the possibility of a layer where he- 
lium mixes with core material. In case (f), we start 
with a mixture of 70/20/10% ^He/i^C/i^O in the re- 
gion 3.98 < r-g < 4.18, i.e. 100 km above and below the 
original core-shell interface. The detonation transcends 
into the core in this case. We surmise that a mixed layer 
on both sides of the interface effectively increases the 
distance from the detonator to the core, thus facilitating 
direct drive. 

In cases (g) and (h), we start with a mixture of 
70/20/10% 4He/i2c/i60 and 50/30/20% ^Re/^^G/^^O, 
respectively, in a 100 km wide region above the original 
interface only. Although the detonation changes when 
it enters the mixed region, becoming hotter and faster, 
it fails to penetrate the core and stay intact in both 
cases. Even with half resolution, the detonation in case 
(g) shows the same instability (unlike case c, for which 
the coarser resolution helped to suppress instabilities in 
the detonation front). A mixture of 80/20% '^He/^^C in 
a 100 km wide region also leads to failure (not shown in 
Fig. 3). 

3.1.1. Heavier Dwarfs 

The above discussion shows that direct drive is not 
likely for Model A if the helium detonation starts where 
the gas is hottest, i.e., the most natural location for a 
detonator. Models B and C both have a denser 1.0 Mq 
core which in principle should detonate more easily than 
the O.8M0 core of Model A. Model C is slightly more 
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Fig. 3. — Density (left, upper colorbar) and temperature (right, lower colorbar) in local 2D simulations of helium detonations hitting 
upon the carbon-oxygen core. The black circles indicate the perimeters of the helium detonation spots at the beginning of each simulation. 
The solid and dashed gray lines represent the mass fraction contours for 49% ^^C and 49% ^''O, respectively. 



compact than Model B. 

Following a similar procedure, we tested for the possi- 
bility of direct ignition with detonators at different alti- 
tudes for these models. The radius of the detonation "hot 
spots" was 20 km and the distance between the core-shell 
interface and the temperature maximum was 149 km and 
106 km in Models B and C, respectively. In both cases, 
direct ignition failed if the detonation was started at this 
altitude. 

If we started the helium detonation 100 km above the 
hottest layer of Model B however, the detonation front 
smoothly transcended into the core. If ignited 50 km 
above the hot layer, the detonation fragmented and died 
upon passing into the core. In Model C, the direct drive 
worked with detonators at altitudes of 100 km, 50 km and 
even 25 km above the hottest layer. At 12 km, however, 
it failed^ These results for direct drive with detonators 
at different altitudes are summarized in Table 2. 

To explore the relevant physics and test the accuracy of 
the CASTRO results, a series of calculations were carried 
out using the KEPLER code. These studies explored 
the critical mass for carbon detonation at the relevant 
densities. A hypothesis we wanted to explore was that 
the critical altitude for a successful propagation of the 
helium detonation into the carbon is given by the radius 
of the critical carbon mass required to generate a self- 
sustaining detonation at the density characterizing the 
carbon-helium interface. As we shall see, this hypothesis 
turned out to be approximately correct. 

In the first study, spheres of 50% carbon, 50% oxy- 
gen were prepared with a constant density of 1.93 x 
lO^gcm"-^, as appropriate to Model B. In a one- 
dimensional spherical region of variable radius, the tem- 

^ Given that this displacement is becoming comparable to the 
size of the detonator itself, it is questionable whether direct drive 
will or will not work for the most massive models. Also the exact 
location of of the hot layer is probably not known to this precision 
due to the one-dimensional mixing length approximation used to 
treat convection. Multi-dimensional simulation of the convection 
during the pre-explosive runaway might provide additional insight. 



perature was given a constant gradient (with respect to 
radius, not mass, even though the code is Lagrangian) 
with values ranging from 2.8 x 10^ K at the center to 
1.0 X 10* K at the edge. This 200 zone region was sur- 
rounded by 300 zones of carbon and oxygen at 10* K. 
All zones had the same radial thickness. The radius of 
the region with the temperature gradient was varied and 
its ability to generate a self-sustaining detonation deter- 
mined. The central value of 2.8 x 10^ K is sufficiently 
high to guarantee carbon burning that was initially su- 
personic. This procedure is very similar to the one pre- 
viously followed by Nienieyer & Woosley (1997), Ropke 
ct al. (2007) and Seitenzahl et al. (2009). 

For a "detonator" radius of 150 km, an initially strong 
detonation decayed away after traversing an additional 
70 km of "cold" carbon and oxygen. For a radius of 200 
km however, the detonation propagated successfully to 
the edge of the grid at 500 km. This is sufficiently far, 
considering that in the real star the density would have 
become higher and the detonation more robust. These 
detonators enclosed masses of 2.73 x 10^* g (150 km) and 
6.47 X lO^^g (200 km) respectively. 

More realistically, the detonator consists of helium, not 
carbon and oxygen and the background is better ap- 
proximated as isobaric rather than constant density. A 
temperature of 2.8 x 10^ K in the carbon also implies 
an unrealistically high energy density not likely to be 
achieved at these low densities. A pure helium sphere 
with variable radius was surrounded by a thick shell of 
50% and 50% ^^O again at 1.93 x io^ gcm'^ with a 
temperature of 1.0 x 10* K. The helium core was com- 
prised of 200 Lagrangian mass shells of equal mass with 
a central temperature of only 1.0 x 10^ K and a temper- 
ature at the edge of 1.0 x 10* K and a constant gradient 
(with respect to radius) in between. The large value of 
central temperature assured supersonic burning and the 
well zoned gradient implied the existence of a region in 
the helium where the phase speed of the burning was 
sonic. The density within the helium core was varied so 
as to provide a constant pressure equal to the pressure 
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in the carbon-oxygen layer. Thus pressure on the entire 
grid was initially 7.25 x 10^^ dyne cm~^. Gravity was ne- 
glected and the pressure was initially held constant by 
applying a boundary pressure equal to the pressure on 
the grid. If the radius of the helium sphere was 211 km 
enclosing 6.37 x 10^* g, the detonation successfully prop- 
agated into the surrounding carbon. Another run with a 
helium sphere radius of 159 km (2.69 x 10^^ g) failed to 
detonate the carbon shell however. 

These results are in good qualitative agreement with 
those obtained using a 2D representation of the full star 
in the CASTRO code (Table 2) . They also agree reason- 
ably well with the pure carbon-oxygen study at constant 
density described above. Apparently the critical altitude 
(or detonator radius) is not very sensitive to the differ- 
ent energy yields of carbon and helium burning. This is 
somewhat surprising, but confirms the hypothesis that 
the critical altitude in the helium layer is one that would 
enclose the critical mass of carbon required for detona- 
tion at the interface density. 

3.2. Shock Collisions and Hot Spots 

To study the sliding helium detonation and the inter- 
nal compressional wave which potentially produces the 
seed for a core detonation, we suppress nuclear burning 
in zones with more than 10% ^^C. This is to prevent 
spurious carbon burning at the core-shell boundary due 
to the mixing of hot helium with cold carbon in a zone, 
and to prevent core detonation before the waves in the 
core have fully converged and produced the highest pos- 
sible temperatures and densities (which we want to de- 
termine) . For the sake of comparability, we here focus on 
the (0.801 + O.143)M0 Model A for the most part. The 
detonation of dwarfs with I.OOM0 cores and lightweight 
helium shells is discussed in §3.2.6. 

3.2.1. One Detonator 

We here consider the case of a single, spherical det- 
onation spot setting off a helium detonation where the 
hottest layer (at rg — 4.27) intersects with the positive 
z-axis in our coordinate system. Runs 1, 4, 5 and 11 all 
represent this case with different numerical resolutions. 
The helium detonation wave wraps around the star (see 
Fig. 4a) and reaches the antipode of the detonator af- 
ter about 1.2 s. For comparison, a sound wave in the 
helium shell would need about 5.6 s to reach the other 
side of the star. This corresponds to a front velocity of 
^ 1.1 X 10^ cms~^. The polar (latitudinal) velocity of the 
ash directly behind the front and the sound speed in the 
ash both are roughly 6.9x 10^ cms~^. Their sum is larger 
than the actual front velocity, consistent with incomplete 
burning rather than a Chapman- Jouguet detonation (cf. 
Sim et al. 2012). The ^^Ni fraction far (-1000 km) be- 
hind the detonation front is about 70%. 

The sliding detonation induces an internal compres- 
sional wave which converges off-center inside the core 
on the z-axis (in general: the axis of symmetry in the 
problem). The average velocity with which the per- 
turbation propagates in negative z-direction from the 
helium ignition point to the final convergence point is 
— 4.0 X lO^cms"^, a little larger than the mean sound 
speed inside the core (3.4 x lO^cms"^). 

The temperatures and densities in some zones of run 1 



are high enough to potentially'^ trigger a core detonation 
where the internal wave finally converges (Niemeycr & 
Woosley 1997; Ropke ct al. 2007; Seitenzahl et al. 2009). 
However, it is likely that the core ignites earlier, at larger 
distance from the center and close to the core-shell in- 
terface on the z-axis, where the internal wave front con- 
verges in the polar direction (but not in radial direction), 
and the converging helium detonation drives a hot, thin 
inflow of ^^Ni in positive z-direction. As a numerical 
experiment, we restarted run 1 without suppressing the 
burning of carbon- rich zones at times t = 1.2 s (when 
the helium detonation converges), 1.3s, 1.4s, and 1.5s 
(when the internal wave converges). In all of these cases, 
a detonation wave in the core forms immediately. 

The values of the density and temperature in the hot 
spot are dependent on the spatial resolution Aa; of the 
grid as well as on the temporal separation At of the con- 
sidered snapshots (for practical reasons, we do not edit 
and analyze the results at every time step), cf. the values 
of temperature and density in the hot spots of runs 1, 4, 
5 and 11 listed in Table 3. However, all results unequiv- 
ocally indicate that the core would detonate. We expect 
that a higher resolution would only increase the likeli- 
hood for detonation. 

3.2.2. Opposite Detonators 

As a limiting case of multi-point ignition, we at first 
consider two detonators at opposite points in the helium 
shell (run 2). While it is unlikely for such detonators 
to be synchronous, this poses a computationally cheap 
2D problem which can be regarded as a numerical ex- 
periment. The internal waves in this case collide on the 
z = plane (the detonators being centered on the z-axis 
at zg = ±4.27), beginning at the outer edge of the core at 
t — 0.60 s and reaching the center at t = 1.05 s. The final 
wave collision in the central part of the star takes place 
on a —1400 km wide disk without generating a discernible 
hot spot. However, temperatures > 10^ K occur on the 
collision plane at large radii (rg — 1.8 ... 3.2). While the 
material is not compressed there, we observe a radially 
inward moving core detonation in a run where the burn- 
ing of carbon is allowed for t > 0.65 s. If carbon burning 
is suppressed, the original core waves are reflected at the 
z = plane. The reflected waves collide on the z-axis 
(symmetrically on both sides of the z = plane), creat- 
ing dense hot spots at zg — ±1.45 that would very likely 
have ignited the core if ignition had failed earlier. 

When one of the detonators (in our case the one on the 
negative side of the z-axis) is delayed by 0.30 s, the two 
helium detonation fronts collide at i = 0.75 s at an angle 
of 112° with respect to the positive z-axis, see Fig. 5. 
As in the synchronous case, the collision of the primary 
internal waves, which here happens on a slightly curved 
disk at z w —600 km, does not generate a hot spot. How- 
ever, a wave reflected toward the negative z-axis pro- 
duces a hot spot that is sufficiently dense to trigger core 

^ Whether a spontaneous detonation initiates depends not only 
on the density, the peak temperature and the size of the region (or, 
equivalently, its mass), but also on the shape of the temperature 
profile (as emphasized by Seitenzahl ct al. 2009). The conditions 
in the hottest zone are listed in Table 3. Since the grid zones 
in our simulations are usually larger than the critical masses for 
detonation, it is not possible to predict with certainty whether the 
conditions in one or a few zones lead to detonation. 
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TABLE 3 

Conditions at the hot spot induced by converging compressional waves in the core 



Run 


Model 




[s] 


[km] 


[s] 


t 

[s] 


z 

[10** cm] 


r° 
[10** cm] 


T 
[10^ K] 


P 

[lO"^ gcm""^] 


Tvnp' 


1 


A 

A 






u.ol 


U.UO 


1 

i.oU 


1 on 
— i.yu 


1 on 
i.yu 


o.zO 


D.Zi 


P 


o 


A 


lOU 




O.Ol 


U.UO 


1 AO 


-1-1 A.^ 
□I i .^0 


1 A^ 


1 0*^ 
1 .yo 


O.OO 


r 


Q 
O 


A 




U.OU 


O.Ol 


U.UO 


1 


1 79 


1 79 


'XA 


A no 

'i.UU 


r 




A 






l^.O 


u.uz 


1 Afi 
± .^o 


1 87 
— i .O 1 


1 9,1 


S 00 

o.uy 


9 fiS 
Z.Do 


P 










O.OD^ 


n 09 
u.uz 




1 09 


1 09 

±.yz 


R1 

y .oi 


9 ^A 

z.o^ 


P 


o 


A 






zo.u 


o^^ 

U.UO 


±.00 


1 Qfl 

— i .yo 


1 Ofi 


00 

o.uy 


A 90 


P 










ZO-U 


nof^ 
u.uuo 


1.0±0 


— i .y i 


1 01 

i .y i 


^ i^i 

O.Ol 


90 

zu.y 


P 


6 


A 


54° 


_ 


25.0 


0.05 


1.40 


-1.76 


1.77 


2.81 


5.66 


P 


7 


A 


54° 


0.15 


25.0 


0.05 


1.45 


-1.74 


1.80 


3.67 


9.51 


P 


8 


A 


90° 




25.0 


0.05 


1.40 


-1.39 


1.54 


2.17 


5.04 


t 


9 


A 


120° 


0.30 


25.0 


0.05 


1.50 


-0.91 


1.96 


2.86 


3.59 


t 


10 


A 


120°, 90°, 100° 


0.30, 0.20 


25.0 


0.05 


1.45 


-0.71 


1.65 


2.19 


5.27 


t 


11'' 


A 






52.1 


0.05 


1.55 


-1.95 


1.95 


2.96 


7.56 


P 


12 


A 


36° 




52.1 


0.05 


1.45 


-1.74 


1.74 


2.19 


4.93 


P 


13 


A 


54° 




52.1 


0.05 


1.45 


-1.74 


1.81 


2.00 


4.24 


P 


14 


E 






6.51 


0.05 


1.40 


-1.70 


1.40 


3.52 


4.54 


P 


15 


E 


180° 




6.51 


0.05 


0.60 





3.48 


3.48 


0.60 


b 


16 


D 






6.51 


0.05 


1.80 


-2.13 


2.13 


2.44 


3.31 


P 


17 


D 


180° 




6.51 


0.05 















^ Angular separation(s) of the detonators in cases with multiple detonators 

Delay of the second (and third) detonator 
° Grid resolution at the highest AMR level that includes the hot spot 
^ Cadence of the considered snapshots 

Distance from the center (spherical radius) 
* Type of listed hot spot: (p) where primary waves converge, (r) where reflected waves collide, (t) where a wave passes through others 
and converges or (b) near the core boundary, where primary waves first collide 
s Restarted with finer zoning in the region where the internal wave converges 
^ Despite this being a 2D problem, this run was done in 3D 



detonation. 

3.2.3. Two Synchronous Detonators 

Before discussing 3D cases of two-point ignition in tlie 
helium, it is worth noting that a single, small aspherical 
detonator is virtually equivalent to a spherical one, re- 
gardless of how strongly deformed it is (as long as it sets 
off a detonation). As the detonation expands with con- 
stant velocity in every direction, the original asphericity 
is quickly lost. To assess the effects of non-axisymmetry 
on the focusing of the internal waves, we therefore con- 
sider cases with two detonators at varying but significant 
separation from one another. 

In cases with two moderately separated, synchronous 
detonators (we consider angular separations of 36°, 54° 
and 90°), the initially separate helium detonation fronts 
collide halfway between the two spots on the geodesic 
connecting them on the core surface, in our case the 
positive z-axis. While the detonations merge, the com- 
bined front proceeds toward the antipode of this collision 
point. The front initially has a head start on the plane 
that connects the two spots with the center of the star 
(the X — plane in our simulations), but the deficit 
becomes smaller as the detonations continue to merge"'. 
As it approaches the antipode, the detonation front has 
the shape of a pointed ellipse. This shape can easily be 
understood considering the combined detonation as a su- 

^ For comparison, consider two merging detonations in a plane, 
starting from the z-axis at ita. With r being the radius of an 
individual detonation, the height of the combined detonation on 
the j/-axis is \/r'^ — ■, asymptotically approaching r (and not r — 
a) at large radii. That is, the two initially separate detonations 
become indistinguishable when their size is much larger than the 
initial separation. 



perposition of two detonations: the pointed ellipse is the 
intersection between two circles (or, more precisely, lines 
of constant latitude, if the ignition points are regarded 
as poles) whose centers are the respective antipodes of 
the ignition points. In a simulation, the detonation fi- 
nally converges when the width of the pointed ellipse 
approaches the size of a grid zone. Theoretically, the 
ratio of the length of the ellipse compared to its width 
increases indefinitely, i.e., the wave always converges in a 
line rather than a point and the difference between hav- 
ing one and two detonators should become larger with 
increasing grid resolution or diffusion length. 

In cases with 36° and 54° detonator separation, the 
region in which the internal waves converge becomes suf- 
ficiently hot and dense to trigger a core detonation, see 
Table 3 and Fig. 4b. We explicitly confirmed this for 
run 13 by switching on carbon burning for t > 1.40 s. 
When two synchronous detonators are 90° apart (run 8), 
there is no hot spot at the site where the two primary 
waves converge at the z-axis (in general; the intersection 
of the two planes of symmetry in this problem). However, 
the waves induced by the two detonators pass through 
one another and form two hot spots at considerable dis- 
tance (660 km) from the z-axis, see Fig. 4c. 

3.2.4. Two Asynchronous Detonators 

We next consider a case of two asynchronous detona- 
tors with an angular separation of 54°, run 7. Note that 
this scenario is less symmetric than the case of two syn- 
chronous detonators, there being only one mirror sym- 
metry in the problem, across the plane which intersects 
both detonators. One of the detonators — the one on the 
+y side in our simulations — is set off 0.15 s (roughly half 
the time needed for the first detonation to reach the sec- 
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Fig. 4. — Snapshots of simulations with one (a) and two synchronous detonators (b: 54° separation, c: 90° separation) during the hehum 
detonation stage (upper panels) and during the presence of a hot spot in the core (lower panels). Panels (b) and (c): the x = plane 
(right half), which contains the two detonators and the center of the star, is symmetric with respect to the y = plane (left half) and vice 
versa; the long ticks at = indicate that these plots show not one but two orthogonal planes. The black contour lines represent, from 
the center outwards: p = 2 X lO^gcm"-^ (long dashed), p = lO^gcm"^ (dotted), Xq = 0.49 (dashed), p = lO^gcm"'' (dot-dashed) and 
p = 10*gcm~^ (dot-dot-dot-dashed). 



blue dashed contour: ''^C^DAS t= 0.750 
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3.00^ 
2.00'^ 
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Fig. 5. — Three snapshots of a 2D simulation with asynchronous detonators on opposite sides (run 3): when the helium detonation 
fronts collide (left), when the internal wave converges (middle) and when reflected waves produce a hot spot (right). The blue dashed line 
indicates the core-shell boundary. Only the inner part of the simulation box is shown. 



ond detonator) after the first one^. Shortly before con- 
verging, the helium detonation front assumes the shape 
defined by the intersection of two circles with different 
radii, see Fig. 6. The same shape can be seen in the in- 
ternal wave front on concentric spherical surfaces. Just 
like in the synchronous case (run 6), a dense hot spot 
forms when the internal wave fronts converge, see Fig. 7 
and the rightmost snapshot in Fig. 6. 

With 120° separation and a delay of 0.30 s, run 9, there 
is no hot spot where the internal waves from the two 
detonators converge. As in the case with synchronous 
detonators at 90° discussed above, the waves from the 
two detonators pass through one another. The wave from 
the first detonator is the first to converge and trigger a 
core detonation. 

3.2.5. Three Asynchronous Detonators 

Unlike the above cases, a setup with three or more 
detonators is in general free of symmetries and requires 

* We implemented this in the code by suppressing the time-step 
update of the zones comprising the detonation spot for said time. 



a 3D simulation of the full star. In run 10, we again 
consider two detonators (A,B) with a separation of 120°, 
one of which (B) is delayed by 0.30 s, exactly as in run 9. 
In addition, a third detonator C, delayed by 0.20 s, is 
placed 90° from detonator A and 100° from detonator B 
on the positive x side. In this case, when the helium det- 
onation converges toward the last patch of unburnt he- 
lium, it roughly has a triangular shape. The same shape 
is present in the internal wave on concentric spherical 
surfaces. There is a distinct hot spot when the inter- 
nal waves converge. Shortly after, a secondary hot spot 
arises as the wave induced by the first detonator con- 
verges again with itself. Fig. 8 shows snapshots of both 
hot spots (right hand panels) and the converging wave 
fronts from which they originate (left hand panels). 

3.2.6. Light Helium Shells 

A detonation in a O.O78M0 helium shell on top of a 
solar mass core (Model E, run 14) is set off by a 25 km 
spherical detonator centered on the hottest helium layer. 
The sliding helium detonation yields a distinct hot spot 
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M-Tg = 1.5 (in core) 

Fig. 6. — Helium shell detonation started by two asynchronous detonators with a 54° separation in the hottest helium layer (run 7). 
The gray, yellow and green isosurfaces indicate the boundary of the CO core, the detonation front and the internal compressional wave, 
respectively. The sketch next to the snapshot at t = 1.05 s (second from right) qualitatively mimics the shape of the converging detonation 
front as the intersection of two circles. Red isosurfaces indicate the hot spot that is generated when the internal wave converges (rightmost 
snapshot). 
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Fig. 7. — Snapshots of the temperature in a simulation with two 
asynchronous detonators at ±27° from the positive 2-axis on the 
X = plane (run 7). The x = plane (right panels) contains the 
detonators and the center of the star, and the orthogonal j/g = 0.46 
plane (left panels) contains a hot spot generated by converging 
internal waves. The detonator on the positive side of the j/-axis 
is delayed by 0.15 s. The black contour lines represent, from the 



10" 



center outwards: p = 2 X 10^ g cm 
(dotted), Xc = 0.49 (dashed), p 
P 

is taken when the internal wave fronts converge 



(long dashed), p = 10 gem 



(dot-dashed) and 

10* gcm~^ (dot-dot-dot-dashed). The snapshot at the bottom 



where the internal wave converges at t 1.40 s. With 
carbon burning turned on from the time when the helium 
detonation converges (t w 1.10 s), the core immediately 
detonates near its surface. This detonation is triggered 
by a thin, hot radial inflow of ^^Ni. As shown in the 
upper panel of Fig. 9, temperatures of 2-4 x 10^ K and 
densities on the order of lO^gcm"'^ are generated along 
this jet. In a simulation with two synchronous, opposite 
detonators, run 15, temperatures > 10^ K are generated 
neither when the wave fronts converge at the center, nor 
along the z-axis when the reflected waves collide (in con- 
trast with Model A, run 2). The only region hot enough 
to light the core in this case is located at the equatorial 
plane near the core's surface, where the collision of the 
internal waves begins. Although there is no hot inflow 




Fig. 8. — Snapshots of the temperature in a simulation with three 
detonators (run 10) on a concentric sphere inside the core which 
contains the hot spot from the three converging waves (upper right 
panel) as well as the hot spot from the wave induced by the first 
detonator alone (lower right panel). 



TABLE 4 

Detonations started in a spherical cap 



altitude'' 

[km] 


arc length 
[km] 


rim density 


rim velocity 
[10* cms-l] 


detonates 





141 


2 X ambient 


8 


yes 





96 


2 X ambient 


8 


no 





233 


ambient 


8 


yes 





141 


ambient 


8 


yes'' 





96 


ambient 


8 


no 





415 


ambient 


1 


yes 





233 


ambient 


1 


no 





415 


ambient 


0.5 


no 





781 


ambient 


0.5 


yes'' 


25 


785 


ambient 


1 


yes 


25 


418 


ambient 


1 


no 



^ Radial distance from the hottest helium layer to the center of 
the cap 

Detonation moves only inward at first and is then reflected back 

comparable to the case with one detonator (note that a 
radial inflow in this case is not a jet), see lower panel in 
Fig 9, the core there ignites in the simulation. 

It is difficult to set off a detonation in a O.O45M0 he- 
lium shell sorrounding a one solar mass core (Model D). 
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Fig. 9. — Conditions near the core boundary as a function of 
distance from the center. Top panel: on the line (negative z-axis) 
where the helium detonation front converges in a run with one deto- 
nator (run 14). Bottom panel; on the helium collision plane (z = 0) 
in a run with two opposite synchronous detonators (run 15). The 
vertical lines mark the initial radius of the core-shell interface. The 
plots represent calculations where carbon burning is suppressed; 
otherwise the core ignites at the earliest time (dot-dot-dot-dashed 
line) in both cases. 
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Fig. 10. — Successful (top) and unsuccessful (bottom) initiation 
of a detonation in the helium shell of Model D with a detonator in 
the form of a spherical cap. The cases shown are the last two in 
Table 4. The color bar ranges from the minimum to the maximum 
temperature in each plot (note the decline of the maximum in case 
of failure). The gray contour line, at 49% ^■^C, marks the core-shell 
interface. 

We did not succeed with a spherical detonator cen- 
tered on the hottest helium layer, whose density is 
7.17 X IG^gcm"^. A detonation does not start even 
with a (perhaps unrealistically) strong initial detonator: 
a central temperature of 3 x 10^ K, linearly decreasing 
to 1.5 X 10^ K at the outer edge, a density three times 
that of the ambient medium and a radius of 25 km (end- 
ing 1 km above the core-shell interface) . Large shock 



velocities (we tried 1.6 x 10^ cms ^ instead of our usual 
8 X 10^ cms^^) at the outer edge of the detonation spot 
did not help either. 

A detonation in the 0.045Mo shell of Model D can be 
set off, however, by a large spherical detonator of radius 
50 km (which is about twice the distance between the 
core-shell interface and the hottest helium layer) cen- 
tered 25 km above the hottest helium layer. The result- 
ing detonation front is more feeble than in all cases with 
heavier helium shells, blowing out at too coarse a reso- 
lution (<2km works, 6.5 km does not), but once started, 
the detonation can propagate around the star. 

We consider the possibility of a thermonuclear runaway 
starting in a (curved) sheet on an equipotential, instead 
of a single point. The size of convective cells in the he- 
lium shell is expected to be on the order of the pressure 
scale height, which is 366 km at the hottest helium layer 
in Model D. The size of the region where the runaway 
starts could be a substantial fraction of this, but is prob- 
ably not larger. As a toy model, we consider detonation 
spots in the form of spherical caps with a thickness of 
50 km. Selected runs are listed in Table 4 and examples 
for a successful and an unsuccessful detonation are shown 
in Fig. 10. As in most of our detonators, the tempera- 
ture decreases linearly from 2 x 10^ K at the center to 
1.8 X 10^ K at the outer edge (on all sides of the cap). 
The outer rims of the caps are half tori, imposed with a 
velocity that increases linearly from the torus center to 
the outer edge. An arc length of 50 km constitutes a lim- 
iting case in which the cap degenerates into a sphere with 
the same properties as our spherical detonators. The cap 
detonator is more powerful than a spherical one. Deto- 
nations can be initiated at ambient density with subsonic 
starter velocities. Smaller sizes require harsher starting 
conditions. Below a size of ^100 km (roughly 1/3 pres- 
sure scale height), the initiation of a detonation appears 
impossible. 

A helium detonation started by a single detonator^, 
run 16, converges after 1.40 s. The corresponding front 
velocity is ~ 9.4 x 10® cms^^ w 4.3 times the sound speed 
in the hottest helium layer. When the internal wave in- 
duced by the helium detonation converges 0.40 s later, 
it forms a distinct hot spot. Without suppression of the 
carbon burning, the core detonates immediately after the 
convergence of the helium detonation, starting at the z- 
axis near the core-shell boundary. In a simulation with 
two synchronous, opposite detonators, run 17, no hot 
spot is formed in the core. With carbon burning turned 
on, the core does not detonate in the simulation. 

3.3. Complete Detonations 

Table 5 lists the yields from helium shell detonations 
in different models, measured at the time when the he- 
lium detonation converges opposite the detonator (we 
only consider single detonators here, not expecting the 
yields from detonations started at multiple points to be 
much different). The lightweight helium shell of Model D 
yields little ^^Ni but relatively large amounts of interme- 
diate mass elements (e.g., *^Ti and "^^Cr). 

^ We note that in order to sustain the helium detonation, the 
shell needs to be more finely resolved than in the other models. We 
here used 1.53 km. A resolution of 6.51km is insufficient, causing 
the helium detonation to blow out quickly. 
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TABLE 5 

Nucleosynthesis yields in the helium shell [Mq] 





Model A 




Model E 




Model D 




12^ 


8.00 X 10- 


4 


5.24 X 10" 


4 


7.02 X 10" 


-4 


"N 


1.11 X IQ- 


10 


2.83 X 10" 


14 


4.24 X 10" 


13 


16o 


3.19 X 10" 


4 


1.93 X 10" 


4 


2.46 X 10" 


-5 






5 


9 79 V 1 n~ 


5 




-5 




9.40 X 10" 


5 


5.30 X 10" 


5 


3.00 X 10" 


-5 




4.46 X 10" 


4 


2.65 X 10" 


4 


2.03 X IQ- 


-4 


32 g 


2.88 X 10" 


4 


1.26 X 10" 


3 


2.17 X 10" 
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Fig. 11. — Mass fractions of elements along the equatorial plane 
(z = 0) in the ejecta of Models E (top panel) and D (bottom 
panel) after complete shell and core detonations, as functions of 
radial velocity (which is approximately proportional to radius) . In 
both panels, only elements with maximum fractions > 10% are 
plotted, respectively. 



We followed the detonation of the core and the expan- 
sion of the ejecta to scales of 10^ km (about 200 times 
the size of the initial dwarf), at which time the expan- 
sion is homologous. The total nickel yields in Models A, 
D and E is O.SSM©, O.52M0 and O.66M0, respectively. 
Mass fractions of selected elements in Models E and D 
are plotted in Fig. 11. The inner parts of the core turn 
into almost pure ^^Ni in both cases, whereas the outer 
parts yield mainly ^®Si and ^^S. The lighter helium shell 
produces large fractions of several intermediate mass el- 
ements instead of ^^Ni, which is the dominant product 



of the heavier shell. 

4. SUMMARY AND DISCUSSION 

We studied the double detonation scenario for Type la 
supernovae by means of 2D and 3D simulations. Starting 
from the results of ID stellar evolution simulations, we 
find that the helium detonation wave is halted at the 
core-shell interface if the detonation starts in the hottest 
layer, which in general is located at some altitude above 
the interface. If the detonation commences at higher 
altitude, the detonation may transcend into the core. For 
direct drive to be successful, the helium detonation must 
involve a critical mass which is sensitive to the density 
at the core-shell interface. A layer of mixed core-shell 
material does not facilitate the direct drive. 

We confirm that the sliding helium detonation induces 
a mildly supersonic compressional wave inside the core 
which, if the detonation is set off at a single point in 
the helium shell, converges to produce an off-center spot 
that is sufficiently hot and dense to light the core as 
well. We tested the robustness of this model for non- 
symmetric initial conditions, considering a series of sce- 
narios in which the helium detonation is started at two 
points and one scenario in which it is started at three 
points. We find that the secondary core detonation is 
quite robust, despite the lack of symmetry. In cases 
where detonators are widely separated (^90°), the hot 
spot is produced by reflected waves or waves that have 
passed through other waves before converging with them- 
selves, rather than the converging primary waves. Only 
in one extreme case — two antipolar, synchronous detona- 
tors in the hehum shefl of a (l.OO2-|-O.O45)M0 dwarf — we 
find no hot spot and no core ignition. 

The generation of a detonation-inducing hot spot is 
robust even in a completely non-symmetric setup with 
three asynchronous detonators. Geometric arguments 
suggest that it is even more robust than a setup with 
two detonators, where the helium detonation and the 
core waves converge in an elongated region. As the elon- 
gation increases with the smallest possible separation of 
the detonation fronts, the difference of having one, two 
or three detonators might be larger at higher resolution. 
Although we have not run cases with more than three 
detonators, we expect that the core ignition would not 
be thwarted by additional detonators. A case with nu- 
merous detonators may be similar to a detonation in a 
spherical shell. 

After the helium detonation has converged, the inter- 
nal waves converge in polar (latitudinal) direction at ever 
decreasing radii, beginning right beneath the core-shell 
interface. In addition, the collision of the helium deto- 
nation front can drive a jet of hot ash into the core. In 
any case, this region has high ignition potential. That 
is, the core may ignite at large radii before the "proper" 
hot spot forms deeper inside as the wave fronts converge 
also in the radial direction. The complete convergence of 
internal waves may only happen in cores that are difficult 
to ignite. 

Our simulations suggest that it is very hard, perhaps 
impossible under realistic assumptions, to initiate a he- 
lium detonation at densities < lO^gcm"'^ (as observed 
in plane-parallel geometry by Townsley ct al. 2012). On 
the other hand, the yields from white dwarfs with thin 
helium shells are most compatible with observed Type la 
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supernovae. Detonations in the helium are more easily 
started in an extended sheet rather than a point. Our 
simulations indicate that the size of such a sheet would 
have to be at least on the order of the size of a convec- 
tive cell in the helium shell. Whether detonation seeds 
of this kind are realistic can ultimately be answered only 
by suitable 3D convection simulations. 
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